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ABSTRACT 

The study of the atmospheres of transiting exoplanets requires a photometric precision, and 
repeatability, of one part in ~10 4 . This is beyond the original calibration plans of current obser¬ 
vatories, hence the necessity to disentangle the instrumental systematics from the astrophysical 
signals in raw datasets. Most methods used in the literature are based on an approximate instru¬ 
ment model. The choice of parameters of the model and their functional forms can some times 
be subjective, causing controversies in the literature. Recently. iMorello et all ( 2014 . 2015 ) have 
developed a non-parametric detrending method that gave coherent and repeatable results when 
applied to Spitzer/IRAC datasets that were debated in the literature. Said method is based on 
Independent Component Analysis (ICA) of individual pixel time-series, hereafter “pixel-ICA’. 
The main purpose of this paper is to investigate the limits and advantages of pixel-ICA on a 
series of simulated datasets with different instrument properties, and a range of jitter timescales 
and shapes, non-stationarity, sudden change points, etc. The performances of pixel-ICA are com¬ 
pared against the ones of other methods, in part icular polyn omial centroid division (PCD), and 
pixel-level decorrelation (PLD) method (Deming et ajJ 20141) . 

We find that in simulated cases pixel-ICA performs as well or better than other methods, and it 
also guarantees a higher degree of objectivity, because of its purely statistical foundation with no 
prior information on the instrument systematics. The results of this paper, together with previ¬ 
ous analyses of Spitzer/IRAC datasets, suggest that photometric precision and repeatability of 
one part in 10 4 can be achieved with current infrared space instruments. 

Subject headings: methods: data analysis - techniques: photometric - planets and satellites: atmospheres 


1. Introduction 

The field of extrasolar planetary transits is one 
of the most productive and innovative subject in 
Astrophysics in the last decade. Transit observa¬ 
tions can be used to me asure the size o f pla nets, 
their orbital parameters ( Seaa gr_fc^fallemOrnelaa 


2003 ), stellar properties dMandel fo Agoll [2002; 


Howartbl l201ll) to study the atmospheres of 


planets ( Brown 2001 : Charbonneau et al. 2002; 


(jBr own 2001J), which is comparable to the effects 
of current instrumental systematics and stellar 


activity ( Berta et al. 20 lit Ballerini et al. 2012 ). 
hence the necessity of testable methods for data 
detrending. In some cases, different assumptions, 
e.g. using different instrumental information or 
functional forms to describe them, leaded to con¬ 
troversial results even from the same datasets^ ex- 


2002 ■ Ago^^O" 200^) , and exomoons dKinninaJ 
MoaaP In particular, the study of planetary at¬ 
mospheres requires a high level of photometric 
precision, i.e. one part in ~10 4 in stellar flux 


amples in the literature are Tinetti et al . 


* — a— i j— - ~ ~ ~'~' i r - t^hrenreichetah _J 2007): loeauneu ei au i|zuuon ; 

Tinetti et al.||2007 ), to detect small planets dMiralda-Escuc Qesert et al.1 ( 2009[_ 2011 ) for the h ot-Jupiter 


Beaulieu et al 


(2007); 

(20081) : 


HD18 9 733b. andlS tevenson et al. d2010l) : Beaulieu et al 


d201lf) ; iKnutson et al.1 1)201 lL 120141) for the warm- 
Neptune GJ436b. Some of these controversies 
are based on Spitzer/IRAC datasets at 3.6 and 
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4.5 fim. The main systematic effect for these 
two channels is an almost regular undulation 
with period ~3000 s, so called pixel-phase ef¬ 
fect, as it is correlated with the relative posi¬ 
tion of the source centroid with respect to a pixel 


center ( Fazio et al.l 120 04; Morales-Cald eron et al 


2006). Conventional parametric techniques cor¬ 
rect for this effect dividing the measured flux by 
a polynomial function of the coordinates of the 
photometric centroid; so me variants may i nclud e 
time-dependen c e (e.g . Stevenson et al.l ( 201Clh : 
Beaulieu et al.1 ( 2011 1). Newer techniques at¬ 
tempt to map the intra-pixel variability at a 
fine-scale level, e.g. adopting spatial weight¬ 
ing f u nctions (Ballard et all l2010t ICowan et al.l 

or interpolating grids 
The results obtained 


2012; Lewis et al 


201 


(Stevenson et al.l 2012U 


with these methods appear to be strongly depen¬ 
dent on a few assumptions, e.g. the degree of 
the polynomial adopted, the photometric tech¬ 
nique, the centroid determination, calibrating in- 


whole observation (e.g. 

Beaulieu et al. 

Diamond-Lowe et al. ( 

2014 

); Zellem et al. 


(120141 11. Also, the very same method, ap¬ 


plied to different observations of the same sys¬ 
tem, often leads to significantly different re¬ 
sults. Non-parametric methods have been pro¬ 
posed to guarantee an higher degree of objec- 


Gibson et al. 2012 

Waldmann 2012: Waldmann et al. 

2013 

: Waldmann 

2014. Morello et al.l (20l4 


20151 ) reanalyzed the 3.6 and 4.5 /im Spitzer/IRAC 


primary transits of HD189733b and GJ436b ob¬ 
tained during the cryogenic regime, so called “cold 
Spitzer” era, adopting a blind source separation 
technique, based on an Independent Component 
Analysis (ICA) of individual pixel time series, 
in this paper called “pixel-ICA”. The results ob¬ 
tained with this method are repeatable over dif¬ 
ferent epochs, and a photometric precision of one 
part in ~10 4 in stellar flux is achieved, with no 
signs of significant stellar vari ability as suggested 


in the previo us _ literature (IDesert et al.l 12011 


Knutson et, al.ll201lt ). The use of ICA to decor 
relate the transit signals from astrophysical and 
instrumental noise, in spectroph otometr ic obser ¬ 
vations, has been propo s ed by Waldmann J2012I) ; 
Waldmann et al.1 (120131) : Waldmann ( 201411 . The 
reason to prefer such blind detrending methods 


is twofold: they require very little, if any, prior 
knowledge of the instrument systematics and as¬ 
trophysical signals, therefore they also ensure a 
higher degree of objectivity compared to meth¬ 
ods based on approximate instrument systematics 
models. As an added value, they give stable re¬ 
sults over several datasets, also in those cases 
where more conventio nal methods have bee n un¬ 
successful. Recently, iDeming et al. ( 2014 ) pro¬ 
posed a different pixel-level decorrelation method 
(PLD) that uses pixel time series to correct for the 
pixel-phase effect, while simultaneously model¬ 
ing the astrophysical signals and possible detector 
sensitivity variability in a parametric way. PLD 
has been applied to some Spitzer/IRAC eclipses 
and synthetic Spitzer data, showing better per¬ 
formances compared to previously published de¬ 
trending methods. 


In this paper, we test the pixel-ICA detrending 
algorithm over simulated datasets, for which in¬ 
strumental systematic effects are fully under con¬ 
trol, to understand better its advantages and lim¬ 
its. In particular, we: 


1. consider some toy-models that can repro¬ 
duce systematic effects similar to those seen 
in Spitzer/IRAC dat£0; 

2. test the pixel-ICA method on simulated 
datasets; 

3. explore the limits of its applicability; 

4. compare its performances with the most 
common parametric method, based on di¬ 
vision by a polynomial function of the cen¬ 
troid (PCD), and the semi-parametric PLD 
method. 


2. Brief outline of pixel-ICA algorithm 


ICA is a statistical technique that transforms a 
set of signals into an equivalent set of maximally 
independent components. It is widely used in a 
lot of different contexts, e.g. Neuroscience, Econo¬ 
metrics, Photography, and Astrophysics, to sepa¬ 
rate the source sig nals present in a set of obser¬ 
vations/recordings ( Hvvarinen et al.1 12001 ). The 


J we do not mean to emulate the Spitzer/IRAC system, but 
rather to study some mechanisms that, in particular con¬ 
figurations, may reproduce in part the effects observed in 
Spitzer/IRAC datasets. 
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underlying assumption is that real signals are lin¬ 
ear mixtures of independent source signals. The 
validity of this assumption for astrophysical ob- 
ser vations has b e en dis cussed with more details 
in iMorello et al.1 ( 2015h (App. A). The major 
strenghts of this approach are: 


1. it requires the minimal amount of prior in¬ 
formation; 

2. even if the assumption is not valid, the 
method is able to detrend in part the source 
signals. 


If additional information is available, some vari¬ 
ants of ICA can be used to obtain better results 


(Iau al et al.ll2002t IStone et alJl2002t iBarriga et al 


201 1). but these are not considered in this paper. 

Pixel-ICA method uses individual pixel time se¬ 
ries from an array to decorrelate the transit signal 
in photometric observations of stars with a tran¬ 
siting planet. The main steps are: 


1. ICA transformation of the time series to get 
the independent components; 


2. identification of the transit com ponent, by 


i nspec tion or other methods (e.g. IWaldmann 

fiH)); 


3. fitting of the non-transit components plus a 
constant on the out-of-transit of the integral 
light-curve (sum over the pixel array); 


4. subtraction of the non-transit components, 
with coefficients determined by the fitting, 
from the integral light-curve; 


5. normalization of the detrended light-curve. 


The norm alized, detrended light -curve is model- 
fitted with Mandel fc Ago! ( 2002 ), which depends 
on several stellar and orbital parameters. We typ¬ 
ically perform a Nelder-Me ad optimisation algo¬ 
rithm (iLagarias et al. 119981) to obtain first esti¬ 


mates of the best parameters of the model, then we 
generate Monte Carlo chains of 20,000 elements to 
sample the posterior distributions (approximately 
gaussians) of the parameters. The updated best 
parameters are the mean values of the chains, and 
the zero-order error bars, cr pa r, 0 ; are their stan¬ 
dard deviations. The zero-order error bars only 


accounts for the scatter in the detrended light- 
curve; they must be increased by a factor that 
includes the uncertainties due to the detrending 
process: 


&nar — 


par, 0 \ 


2 

a ICA 


(1) 


where a par is the final parameter error bar, ao is 
the square root of the likelihood’s variance (ap¬ 
proximately equal to the standard deviation of 
residuals), and ar c a is a term as s ociated to th e de¬ 
trending process. Morello et al. ( 2014 . 2015h sug¬ 
gest the following formula for ajcA ■ 


2 _ f 2 
a ICA ~ J 


cdlSR, +< 


( 2 ) 


where ISR is the so-called Interference-to-Signal- 
Ratio matrix, Oj are the coefficients of the non¬ 
transit-components, to is their number, cr n t c -fit 
is the standard deviation of residuals from the 
theoretical raw light-curve, out of the transit, / 
is the normalising factor for the detrended light- 
curve. The sum on the left takes into account 
the precision of the components extracted by the 
algorithm; a n tc-fit indicates how well the linear 
combination of components approximates the out- 
of-transit. The MULTICOMBI code, i.e. the al¬ 
gorithm that we use for the ICA transformation, 
provides two Interference-to-Signal-Ratio matri¬ 
ces, ISR^ and ISR." ' 4 5 , associated to the sub¬ 
algorithms EFICA and WASOBI, respectively. 
Two approaches has been suggested to derive a 
single Interference-to-Signal-Ratio matrix: 


ISR = 


ISR ef + ISR 1- " 4 
2 


(3) 


ISRij = min (iSH/'/'.lSR! 1 ; 1 ) (4) 

Eq. [3] is a worst-case estimate, while Eq. []] takes 
into account the outperforming separation capa¬ 
bilities of MULTICOMBI compared to EFICA 
and WASOBI. We adopt Eq. [I] throughout this 
paper, but results obtained with both options are 
reported in Tab. [7] [31 0 and EH In most cases 
the differences are negligible. 
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3. Instrument simulations 

3.1. Instrument jitter only 

We consider an ideal transit light-curve with 
parameters reported in Tab. [l] sampled at 8.4 
s over 4| hr, totaling 2001 data points, symmet¬ 
ric with respect to the transit minimum. Fig. [L] 


pixel. The same is valid for the flux integrated over 
a small aperture, compared to the width of the 
PSF. Also note that short time scale fluctuations 
are present as a sampling effect. Neither of these 
variations are observed with a large aperture that 
includes all the photons at any time. 

3.2. Source Poisson noise 


Table 1: Transit parameter values adopted in all 
simulations: p = r p /R s is the ratio of planetary 
to stellar radii, ao = a/R s is the orbital semima¬ 
jor axis in units of stellar radius, i is the orbital 
inclination, e is the eccentricity, P is the orbital 
period, 71 and 7 ? are quadratic limb darkening 
coefficients (IHowarthl ,2011). 


p 

0.15500 

ao 

9.0 

i 

85.80 

e 

0.0 

p 

2.218573 days 

71,72 

0.0, 0.0 


shows the ideal light-curve. To each data point 
we associate a number of photons proportional to 
the expected flux, in particular 50,000 photons in 
the out-of-transit. We generate random gaussian 
coordinates for each photon, representing their po¬ 
sitions on the plane of the CCD. Finally, we add 
a grid on this plane: each square of the grid rep¬ 
resents a pixel, and the number of photons into a 
square at a time is the read of an individual pixel 
in absence of pixel systematics. To simulate the 
effect of instrumental jitter, we move the gridlines 
from one data point to another (it is equivalent to 
shift the coordinates of the photons). 



In a real case scenario, the detected photons 
from the astropliysical source are distributed ac¬ 
cording to a poissonian function with fluctuations 
proportional to the square root of the expected 
number of photons. The value of the proportion¬ 
ality factor is specific to the instrument. Poisson 
noise in the source is a natural limitation to the 
photometric precision that can be achieved for any 
astropliysical target. To include this effect in sim¬ 
ulations, we added a random time series to the as- 
trophysical model, multiplied by different scaling 
factors, then generated frames with total number 
of photons determined by those noisy models. 


3.3. The effect of pixel systematics 


Spitzer/IRAC datasets for channels 1 and 2 
show temporal flux variations correlated to the 
centroid posit ion, independently from the aperture 
select ed (Fazio et al .2Q0ll Morales-Calderon et al. 


20061 ). Here we consider two effects that, com¬ 


bined to the instrument jitter, can produce this 
phenomenon: 


1. inter-pixel quantum efficiency variations, 
simulated multiplying the photons in a pixel 
by a coefficient to get the read, being the 
coefficients not identical for all the pixels; 

2. intra-pixel sensitivity variations, simulated 
by assigning individual coefficients depen¬ 
dent on the position of the photon into the 
pixel. 

3.4. Description of simulations 


Fig. 1. - Referent transit light-curve adopted in 
simulations. 


We performed simulations for two values of 
(gaussian) PSF widths, <JpsF'- 


Fig. [2] shows the jitter effect at different levels, 
i.e. individual pixels, and small and large clusters. 
Pixel flux variations are related to the changing 
position of the PSF: the flux is higher when the 
centre of the PSF is closer to the centre of the 


• &psf =1 p.u. (pixel side units); 

• &psf =0.2 p.u. 
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Tlme ( s ) xio 4 


Fig. 2.— Top panel: Representation of one simulated frame onto the focal plane. Bottom panels: simulated 
time series associated to selected individual pixels (blue and red), a centered 5x5 array (green), and a centered 
9x9 array (black). The centroid is assumed oscillating in the direction indicated by the double-headed arrow, 
with a sinusoidal pattern (sinl, see Tab. [2]). 
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The two sets of frames differ only in the scaling 
factor in photon coordinates, therefore no rela¬ 
tive differences are attributable to random gen¬ 
eration processes. As a comparison, the nominal 
PSF widths for Spit zer/ IRAC channels are in the 


range 0.6-0.8 p.u. (IFazio et al. 20041) . By com¬ 


paring two more extreme cases, we observe the 
impact of the PSF width on pixel systematics and 
their detrending. We also tested eleven Poisson 
noise amplitudes, linearly spaced between 0 and 
^50,000. 

We simulated several jitter time series as detailed 
in Tab. [2]and Fig. [3] For each case, we adopted: 


1. a random generated quantum efficiency map 
with standard deviations of ~10 -2 , to sim¬ 
ulate the inter-pixel variations; 

2. a non-uniform response function for each 
pixel, i.e. 1—0.Id, where d is the distance 
from the centre of the pixel, to simulate the 
intra-pixel variations. 


These values are of the same order as the nom¬ 
inal pixel-to-pixel accuracy of flat-fielding and 
sub-pixel response fo r Sp itzer/IRAC channels 1 
and 2 ( Fazio et al.112004! ) , but slightly larger to 
better visualize their effects. Finally, we add 
white noise time series at an arbitrary level of 5 
counts/pixel/data point for most cases, to sim¬ 
ulate high-frequency pixel noise (HFPN). It is 
worth to note that the same quantum efficiency 
maps and noise time series have been adopted for 
all the simulations with pixel arrays of the same 
size, to minimize possible aleatory effects when 
comparing the results. 
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Table 2: List of jitter time series adopted in simulations and their properties. 


Abbr. 

Shape 

Peak-to-peak amplitude (p.u.) 

Period (s) 

sinl 

sinusoidal 

0.6 

4014.6 

cost 

cosinusoidal 

0.6 

4014.6 

sin2 

sinusoidal 

0.6 

3011.0 

cos2 

cosinusoidal 

0.6 

3011.0 

sin3 

sinusoidal 

0.6 

2007.3 

cos3 

cosinusoidal 

0.6 

2007.3 

sawl 

~sawtooth 

0.6 

2990.4 

sawlvl 

~sawtooth 

variable 

2990.4 

sawlv2 

~sawtooth 

variable 

2990.4 

sawlv3 

~sawtooth 

decreasing 

2990.4 

sawlvfl 

~sawtooth 

0.6 

variable 

sawlvf2 

~sawtooth 

0.6 

decreasing 

jump04c 

Heaviside step 

0.4 

mid-transit discontinuity 



-1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 

Time (s) w 



Time (s) 


Fig. 3. Top panel: jitter time series sawl (black), sawlvl (ecru, cross markers), sawlv2 (green, triangles), 
and sawlv3 (cyan, empty circles); markers are represented every 20 data points for reasons of visibility. 
Bottom panel: jitter time series sawlvfl (orange), and sawlvf2 (grey). The other jitter time series are not 
reported, since their representations are obvious (see Tab. 0. 
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4. Results 

After having generated the simulated raw 
datasets, we applied both pixcl-ICA and PCD 
detrending techniques to evaluate their relia¬ 
bility and robustness in those contexts. More 
specifically, the polynomial decorrelating function 
adopted is: 

p(x,y) = a 1 +a 2 (x-x)+a 3 (y-y)+a 4 (x-x) 2 +a 5 {y-y) 2 

(5) 

where x and y are the centroid coordinates, calcu¬ 
lated as weighted means of the counts per pixel, 
x and y are the mean centroid coordinates on the 
out-of-transit, and cq coefficients are fitted on the 
out-of-transit. In several cases we tested some 
variants of PCD, i.e. higher order polynomials, 
cross terms, etc. 

Since the main purpose of this paper is to test 
the ability of pixel-ICA method to detrend instru¬ 
ment systematics, we first report the detailed re¬ 
sults for all the configurations obtained in the limit 
of zero stellar noise. We show in particular (Sec. 
rnirra and roii : 

• the simulated raw light-curves and the cor¬ 
responding detrended ones; 

• the root mean square (rms) from the theo¬ 
retical transit light-curve, before and after 
the detrending processes; 

• the residual systematics in the detrended 
light-curves; 

• the planetary, orbital, and stellar parameters 
estimated by fitting the light-curves; 

• for a subsample of cases, the results of full 
parameter retrieval, including error bars. 

The impact of stellar poissonian noise at different 
levels is studied for a few representative cases (Sec. 
14.411 . A variant of pixel-ICA algorithm, and PLD 
detrending method are discussed in Sec. 14.61 and 

m 

4.1. Case I: inter-pixel effects, large PSF 

Fig. [4] shows the raw light-curves simulated 
with apsF =1 p.u., and inter-pixel quantum ef¬ 
ficiency variations over 9x9 array of pixels, and 
the correspondent detrended light-curves obtained 


with pixel-ICA and PCD methods. This array 
is large enough that the observed modulations 
are only due to the pixel effects (see Sec. 14.11) . 
Tab. [3] reports the discrepancies between the de¬ 
trended light-curves and the theoretical model. 
The amplitude of residuals after PCD detrending, 
i.e. 3.0 xlO -4 for the selected binning, equals the 
HFPN level (see Sec. 13.411 . while for pixel-ICA 
detrended data it is smaller by a factor ~l/3. 
We note that two independent components are 
enough to correct for the main instrument sys¬ 
tematics within the HFPN level, the other compo¬ 
nents slightly correct for the residual pixel system¬ 
atics. Interestingly, a similar behaviour h as been 


obser ved for real Spitzer datasets (IMorello et al 


[mi, but with the best fit model instead of the 
(unknown) theoretical one. Fig. [5] shows how 
the residuals scale for binning over n points, with 
1 < n < 10. This behaviour suggests a high level 
of temporal structure in raw data, which is not 
present in ICA-detrended light-curves. Some sys¬ 
tematics are still detected in residuals obtained 
with the parametric method. We checked that 
the empirical centroid coordinates are accurate 
within 0.006 p.u. on average, and higher order 
polynomial corrections lead to identical results. 
Fig. [6] shows the transit parameters retrieved from 
the detrended light-curves; in a few representative 
cases, we calculated the error bars as detailed in 


Morello et al.1 (12015 ), and Sec. [2] in this paper. 


Numerical results are reported in Tab. |T| 







Raw 9x9, sinl, inter, PSF = 1 


Raw 9x9, cosl, inter, PSF = 1 




Raw 9x9, cos2, inter, PSF = 1 



Raw 9x9, cos3, inter, PSF = 1 



Raw 9x9, sawlvl, inter, PSF = 1 



Raw 9x9, saw1v3, inter, PSF = 1 



Raw 9x9, saw1vf2, inter, PSF = 



Raw 9x9, sinl, inter, PSF = 1 


0.99 

0.98 


“Sr 


0.99 

0.98 


“ST 


0.99 

0.98 


°#05 


0.99 

0.98 




0.99 

0.98 




0.99 

0.98 




I PCD] 

I ICA I 


| PCD] 

I ' 1CA I 


| PCD] 

I 'CA | 


| PCD] 

I 'CA | 


| PCD] 

I 'CA | 


| PCD] 

I ICA | 



Fig. 4.— Left panels: (blue) raw light-curves simulated with upsf = lp.w., and inter-pixel quantum efficiency 
variations over 9x9 array of pixels. Right panels: detrended transit light-curves obtained with (green ‘x’) 
polynomial centroid fitting method, and (red dots) pixel-ICA method. All the light-curves are binned over 
10 points, except those in the bottom right, to make clearer visualization of the systematic effects. 
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Table 3: Root mean square of residuals between the light-curves and the theoretical model for simulations 
with crpsF = 1 p.u., and inter-pixel quantum efficiency variations over 9x9 array of pixels; in particular they 
are calculated for the raw light-curves, light-curves detrended with pixel-ICA, and PCD method, binned over 
10 points. 


Jitter 

rms (raw — theoretical) 

rms (ICA — theoretical) 

rms (PCD — theoretical) 

sinl 

6.5xl0 -4 

2.0xl0 -4 

3.0xl0~ 4 

cosl 

6.9xl0 -4 

2.1X10" 4 

3.2x 10 -4 

sin2 

6.7xl0 -4 

2.0xl0 -4 

3.Ox 10 -4 

cos 2 

6.5xl0" 4 

2.1xl0" 4 

3.1xl0" 4 

sin3 

6.6xl0 -4 

2.0xl0" 4 

3.Ox 10 -4 

cos3 

6.6xl0 -4 

2.0xl0" 4 

3.1xl0" 4 

sawl 

5.4xl0 -4 

2.1xl0" 4 

3.1xl0" 4 

sawlvl 

6.0xl0 -4 

2.1xl0" 4 

3.Ox 10 -4 

sawlv2 

5.6xl0 -4 

2.1xl0" 4 

3.1xl0" 4 

sawlv3 

5.9xl0 -4 

2.1xl0 -4 

3.1x 10 -4 

sawlvfl 

5.5xl0" 4 

2.1xl0 -4 

3.Ox 10 -4 

sawlvf2 

5.3xl0 -4 

2.0xl0 -4 

3.Ox 10 -4 

jump04c 

6.9xl0" 4 

1.9xl0 -4 

3.Ox 10 -4 


9x9 array, inter, PSF = 1 



i i i i i i i i 


1 


! i , 


« ‘'ii. 


‘ • ! i ■ 



---pure white noise 


• raw - theoretical 


* ICA - theoretical 

* * * , S 

PCD - theoretical 

,,,, , , 


1 


2 3 4 56789 10 

binning 


Fig. 5.— Root mean square of residuals for light-curves binned over 1 to 10 points, scaled to their non-binned 
values. The simulations were obtained with <tpsf =1 p.u., 9x9 array, and inter-pixel effects. The dashed 
black line indicates the expected trend for white residuals, blue dots are for normalized raw light-curves, red 
V are for pixel-ICA detrendend light-curves, and green ‘x’ for PCD detrended light-curves. 
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0.1560 


0.1555 


Q-0.1550 


0.1545 


0.1540 


9.5 


8.5 


87.0 


86.5 


- 86.0 


85.5 


85.0 


X 

• ICA 
x PCD 
-■■ref. value 


i 2 3 4 5 6 7 8 9 10 11 12 13 

Jitter cases 



• ICA 
x PCD 
--■ref. value 


i 2 3 4 5 6 7 8 9 10 11 12 13 

Jitter cases 



• ICA 
x PCD 
-■■ref. value 


i 2 3 4 5 6 7 8 9 10 11 12 13 



Jitter cases 


Fig. 6.- Top panel: best estimates of the planet-to-star radii ratio, p = r p /R s , for detrended light-curves 
with (red dots) pixel-ICA, and (green x) PCD method (apsF =1 p.u., inter-pixel effects over 9x9 array). 
Error bars are reported for a few representative cases of jitter signal, i.e. sinl, cosl (chosen as examples of 
periodic functions with different phasing), sawlv3 (example with non-stationary amplitude), sawlvf2 (non¬ 
stationary frequency), and jump04c (sudden change). Middle panel: the same for the orbital semimajor axis 
in units of the stellar radius, ao = a/R s . Bottom panel: the same for the orbital inclination, i. 
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For the same configuration, i.e. <tpsf =1 
p.u., inter-pixel effects, we investigated the con¬ 
sequences of considering a smaller array (5x5), 
which does not include the whole PSF. Fig. Q] 
shows the raw light-curves, and the correspondent 
detrended ones, obtained with the two methods. 
Tab. |4] reports the discrepancies between those 
light-curves and the theoretical model. The dis¬ 
crepancies are higher than for the larger pixel- 
array by a factor > 2 (for the raw light-curves), 
because of the additional effect. After pixel-ICA 
detrending, the discrepancies are reduced by a fac¬ 
tor ~5 (for the selected binning) in most cases, and 
~13 for the ‘jump04c’, while the performances of 
the parametric method are case dependent, and 
discrepancies are reduced by a factor between 2 
and 7 in most cases, and also ~13 for ‘jump04c’. 
In all cases, the final discrepancies are higher than 
the HFPN level, i.e. ~1.7x 10 -4 for the 5x5 array. 
Fig. [8] shows how the residuals scale for binning 
over n points, with 1 < n < 10. The temporal 
structure due to jitter effect is dominant in raw 
data, but little traces of this behaviour are present 
after pixel-ICA detrending. Even for this aspect, 
the performances of the parametric method are 
case dependent. Fig. [2] shows the transit param¬ 
eters retrieved from detrended light-curves; in a 
few representative cases, we calculated the error 
bars. Numerical results are reported in Tab. [5] In 
conclusion, the choice of a non-optinral pixel array 
introduces additional systematics, that worsen the 
parameter retrieval, but it is quite remarkable that 
the pixel-ICA technique gives consistent results in 
most cases, whereas the parametric technique ap¬ 
pears to be less robust. 
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Raw 5x5, sinl, inter, PSF = 1 Detrended Raw 5x5, cosl, inter, PSF = 1 Detrended 



Fig. 7. Left panels: (blue) raw light-curves simulated with crpsF = 1 p.u., and inter-pixel quantum efficiency 
variations over 5x5 array of pixels. Right panels: detrended transit light-curves obtained with (green ‘x’) 
polynomial centroid fitting method, and (red dots) pixel-ICA method. All the light-curves are binned over 
10 points to make clearer visualization of the systematic effects. 
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Table 4: Root mean square of residuals between the light-curves and the theoretical model for simulations 
with crpsF = 1 p.u., and inter-pixel quantum efficiency variations over 5x5 array of pixels; in particular they 
are calculated for the raw light-curves, light-curves detrended with pixel-ICA, and PCD method, binned over 
10 points. 


Jitter 

rms (raw — theoretical) 

rms (ICA — theoretical) 

rms (PCD — theoretical) 

sinl 

1.5xl0 -3 

3.1xl0 -4 

2.6xl0“ 4 

cosl 

1.4xl0 -3 

3.4xl0 -4 

8.1xl0" 4 

sin2 

1.5xl0 -3 

3.2xl0 -4 

2.8xl0" 4 

cos 2 

1.5xl0" 3 

3.2xl0 -4 

6.2xl0" 4 

sin3 

1.5xl0 -3 

3.1xl0" 4 

2.7xl0" 4 

cos3 

1.4xl0 -3 

3.2xl0" 4 

4.9x 10 -4 

sawl 

1.7xl0 -3 

3.1xl0" 4 

3.5x 10 -4 

sawlvl 

1.7xl0 -3 

2.9xl0" 4 

2.5xl0" 4 

sawlv2 

1.6xl0 -3 

3.6xl0" 4 

3.7x 10 -4 

sawlv3 

1.7xl0 -3 

3.1xl0 -4 

4.Ox 10 -4 

sawlvfl 

1.7xl0" 3 

3.1xl0 -4 

2.9xl0" 4 

sawlvf2 

1.6xl0 -3 

3.2xl0 -4 

2.6xl0" 4 

jump04c 

3.3xl0 -3 

2.6xl0 -4 

2.6x 10 -4 


5x5 array, inter, PSF = 1 


i 


--pure white noise 

• raw - theoretical 

* ICA - theoretical 
PCD - theoretical 



2 


3 4 5 6789 10 

binning 


Fig. 8.— Root mean square of residuals for light-curves binned over 1 to 10 points, scaled to their non-binned 
values. The simulations were obtained with <tpsf =1 P-u., 5x5 array, and inter-pixel effects. The dashed 
black line indicates the expected trend for white residuals, blue dots are for normalized raw light-curves, red 
V are for pixel-ICA detrendend light-curves, and green ‘x’ for PCD detrended light-curves. 
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Fig. 9.— Top panel: best estimates of the planet-to-star radii ratio, p = r p /R s , for detrended light-curves 
with (red dots) pixel-ICA, and (green ‘x’) PCD method (ctpsf =1 p.u., inter-pixel effects over 5x5 array). 
Error bars are reported for representative cases of jitter signal, i.e. sinl, cosl, sawlv3, sawlvf2, and jump04c. 
Middle panel: the same for the orbital semimajor axis in units of the stellar radius, ao = a/R s . Bottom 
panel: the same for the orbital inclination, i. 


15 






















4.2. Case II: inter-pixel effects, narrow 
PSF 

Fig. [TO] shows the raw light-curves simulated 
with (JpsF =0.2 p.u., and inter-pixel quantum effi¬ 
ciency variations over 5x5 array of pixels, and the 
correspondent detrended light-curves, obtained 
with the two methods considered in this paper. 
The array is large enough that observed modu¬ 
lations are due only to the pixel effects. Tab. 
[5] reports the discrepancies between those light- 
curves and the theoretical model. Again, the dis¬ 
crepancies after pixel-ICA detrending are below 
the HFPN level, i.e. 1.5xl0~ 4 for the selected 
binning, outperforming the parametric method 
by a factor 2-3. We also noted that the em¬ 
pirical centroid coordinates may differ from the 
“true centroid coordinates” up to 0.15-0.33 p.u., 
which is an error comparable with the jitter am¬ 
plitude. This is not surprising, given that the PSF 
is undersampled, being much narrower than the 
pixel size. Despite the large errors in empirical 
centroid coordinates, in some cases, the discrep¬ 
ancies between PCD detrended light-curves and 
the theoretical model are at the HFPN level, and 
slightly larger in other cases. Fig. ITT! shows how 
the residuals scale for binning over n points, with 
1 < n < 10. A significant temporal structure is 
present in the raw data, but not in the pixel-ICA 
detrended light-curves, while the performances of 
the parametric method are case dependent. Fig. 
M shows the transit parameters retrieved from 
detrended light-curves; in representative cases, we 
calculated the error bars. Numerical results are 
reported in Tab. El 
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Raw 5x5, sinl, inter, PSF = 0.2 Detrended Raw 5x5, cosl, inter, PSF = 0.2 Detrended 



Fig. 10. - Left panels: (blue) raw light-curves simulated with crpsF =0.2 p.u., and inter-pixel quantum 
efficiency variations over 5x5 array of pixels. Right panels: detrended transit light-curves obtained with 
(green ‘x’) polynomial centroid fitting method, and (red dots) pixel-ICA method. All the light-curves are 
binned over 10 points to make clearer visualization of the systematic effects. 
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Table 5: Root mean square of residuals between the light-curves and the theoretical model for simulations 
with apsF = 0.2 p.u., and inter-pixel quantum efficiency variations over 5x5 array of pixels; in particular 
they are calculated for the raw light-curves, light-curves detrended with pixel-ICA, and PCD method, binned 
over 10 points. 


Jitter 

rms (raw — theoretical) 

rms (ICA — theoretical) 

rms (PCD — theoretical) 

sinl 

7.2xl0 -4 

7.2xl0 -5 

1.9xl0 -4 

cosl 

7.0xl0" 4 

7.7xl0" 5 

2.8xl0" 4 

sin2 

7.1xl0" 4 

7.4xl0 -5 

1.8x 10 -4 

cos 2 

7.4xl0" 4 

7.4xl0 -5 

2.4x 10 -4 

sin3 

7.1xl0 -4 

7.8xl0 -5 

1.8x 10 -4 

cos3 

7.0xl0 -4 

7.5xl0 -5 

2.Ox 10 -4 

sawl 

2.6xl0 -4 

6.7xl0" 5 

1.5xl0" 4 

sawlvl 

2.3xl0 -4 

6.8xl0 -5 

1.5xl0" 4 

sawlv2 

2.3xl0 -4 

6.6xl0 -5 

1.5xl0" 4 

sawlv3 

2.3xl0 -4 

6.8xl0 -5 

1.5xl0" 4 

sawlvfl 

2.3xl0" 4 

6.6xl0 -5 

1.5xl0" 4 

sawlvf2 

2.4x 10 -4 

6.8xl0 -5 

1.5xl0" 4 

jump04c 

7.6xl0" 4 

7.4xl0 -5 

1.5xl0" 4 


5x5 array, inter, PSF = 0.2 


"Nk 


-pure white noise 
raw - theoretical 
ICA - theoretical 
PCD - theoretical 


I 


'•-I.. 
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binning 
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Fig. 11. Root mean square of residuals for light-curves binned over 1 to 10 points, scaled to their non- 
binned values. The simulations were obtained with upsf =0.2 p.u., 5x5 array, and inter-pixel effects. 
The dashed black line indicates the expected trend for white residuals, blue dots are for normalized raw 
light-curves, red V are for pixel-ICA detrendend light-curves, and green ! x’ for PCD detrended light-curves. 
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Fig. 12.— Top panel: best estimates of the planet-to-star radii ratio, p = r p /R Sl for detrended light-curves 
with (red dots) pixel-ICA, and (green ‘x’) PCD method (<Jpsf =0.2 p.u., inter-pixel effects over 5x5 array). 
Error bars are reported for representative cases of jitter signal, i.e. sinl, cosl, sawlv3, sawlvf2, and jump04c. 
Middle panel: the same for the orbital semimajor axis in units of the stellar radius, a o = a/R s . Bottom 
panel: the same for the orbital inclination, i. 
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4.3. Case III: intra-pixel effects 

For simulations with apsF =1 p.u., the effect of 
intra-pixel sensitivity variations is negligible, i.e. 
~ 10 ^ 5 , unless we consider unphysical or very un¬ 
likely cases, where the quantum efficiency can as¬ 
sume both positive and negative values in a pixel, 
or it is zero for a significant fraction of the area 
of the pixel. Intra-pixel effects become significant 
when the PSF is narrower, therefore we analyzed 
only the relevant simulations with apsF =0.2 p.u. 

Fig. [CJ] shows the raw light-curves simulated 
with crpsF =0.2 p.u., and intra-pixel quantum ef¬ 
ficiency variations over 5x5 array of pixels, and 
the correspondent detrended light-curves. The ar¬ 
ray is large enough that the observed modulations 
are only due to the pixel effects. Tab. [G] reports 
the discrepancies between those light-curves and 
the theoretical model. The pixel-ICA technique 
reduces the dicrepancies by a factor 4-8 (for the 
selected binning) for the first 12 jitter series, and 
by a factor 83 for the case ‘jump04c’, outperform¬ 
ing the parametric method by a factor 2-4. Fig. 
m shows how the residuals scale for binning over 
n points, with 1 < n < 10. In this case, the 
temporal structure is preserved in all detrended 
light-curves, except for the case ‘jump04c’, which 
means that both methods have some troubles to 
decorrelate intra-pixel effects. Fig. [15] shows the 
transit parameters retrieved from detrended light- 
curves; in representative cases, we calculated the 
error bars. Detailed numerical results are reported 
in Tab. [TO] While in some cases the paramet¬ 
ric method may perform better than pixel-ICA, if 
adopting higher order polynomials, in some other 
cases higher order polynomials lead to worse re¬ 
sults than lower order polynomials. Higher order 
polynomials might approximate better the out- 
of-transit baseline, but the residuals to the the¬ 
oretical light-curve are not necessarily improved, 
hence the correction can bias the astrophysical sig¬ 
nal. The pixel-ICA method is less case depen¬ 
dent. We do not suggest that the same conclusions 
are necessarily valid for Spitzer, since our simu¬ 
lations have stronger systematics and much nar¬ 
rower PSFs. It is quite remarkable that, though 
the systematics are not well decorrelated, the pa¬ 
rameter retrieval gives the correct results within 
the error bars. 
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Raw 5x5, sinl, intra, PSF = 0.2 Detrended Raw 5x5, cosl, intra, PSF = 0.2 Detrended 



Fig. 13.- Left panels: (blue) raw liglit-curves simulated with crpsF =0.2 p.u., and intra-pixel quantum 
efficiency variations over 5x5 array of pixels. Right panels: detrended transit light-curves obtained with 
(green ‘x’) polynomial centroid fitting method, and (red dots) pixel-ICA method. All the light-curves are 
binned over 10 points to make clearer visualization of the systematic effects. 
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Table 6: Root mean square of residuals between the light-curves and the theoretical model for simulations 
with crpsF = 0.2 p.u., and intra-pixel quantum efficiency variations over 5x5 array of pixels; in particular 
they are calculated for the raw light-curves, light-curves detrended with pixel-ICA, and PCD method, binned 
over 10 points. 


Jitter 

rms (raw — theoretical) 

rms (ICA — theoretical) 

rms (PCD — theoretical) 

sinl 

3.5xl0 -3 

5.1xl0 -4 

1.2xl0" 3 

cosl 

3.5xl0" 3 

5.1xl0 -4 

1.8xl0" 3 

sin2 

3.4xl0" 3 

4.9xl0 -4 

1.2xl0" 3 

cos 2 

3.5xl0 -3 

5.0xl0 -4 

1.5xl0" 3 

sin3 

3.4xl0 -3 

5.0xl0 -4 

1.2x 10 -3 

cos3 

3.4xl0 -3 

5.0xl0" 4 

1.4xl0" 3 

sawl 

3.4xl0 -3 

7.4xl0" 4 

1.8xl0" 3 

sawlvl 

3.8xl0 -3 

7.3xl0" 4 

1.7xl0" 3 

sawlv2 

3.6xl0 -3 

6.8xl0 -4 

1.7xl0" 3 

sawlv3 

3.7xl0 -3 

7.3xl0" 4 

1.7xl0" 3 

sawlvfl 

3.2xl0" 3 

6.6xl0 -4 

1.7xl0" 3 

sawlvf2 

3.0xl0" 3 

6.1xl0" 4 

1.5xl0" 3 

jump04c 

6.8xl0 -3 

8.2xl0 -5 

1.6xl0" 4 


5x5 array, intra, PSF = 0.2 
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Fig. 14.— Root mean square of residuals for light-curves binned over 1 to 10 points, scaled to their non- 
binned values. The simulations were obtained with apsF =0.2 p.u., 5x5 array, and intra-pixel effects. 
The dashed black line indicates the expected trend for white residuals, blue dots are for normalized raw 
light-curves, red V are for pixel-ICA detrendend light-curves, and green ! x’ for PCD detrended light-curves. 
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Fig. 15.— Top panel: best estimates of the planet-to-star radii ratio, p = r p /R s , for detrended light-curves 
with (red dots) pixel-ICA, and (green l x’) PCD method (apsF =0.2 p.u., intra-pixel effects over 5x5 array). 
Error bars are reported for representative cases of jitter signal, i.e. sinl, cosl, sawlv3, sawlvf2, and jump04c. 
Middle panel: the same for the orbital semimajor axis in units of the stellar radius, ao = a/R s . Bottom 
panel: the same for the orbital inclination, i. 
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4.4. The effect of astrophysical source 
Poisson noise 

We tested the effect of astrophysical source 
Poisson noise at different levels. It is well ex¬ 
plained by the results obtained with: 

1 . a psf =1 p.u., inter-pixel variations over 
9x9 array, jitter ‘sinl 

2. &psf =0.2 p.u., intra-pixel variations over 
5x5 array, jitter ‘sinl’ and ‘cosl’; 

Fig. [16] shows the raw and detrended light-curves 
obtained for the listed cases with two finite val¬ 
ues of Signal-to-Noise Ratio, SNR = 447 (inter¬ 
mediate) and SNR = 224 (lowest). Fig. [TT] shows 
how the residuals scale for binning over n points, 
with 1 < n < 10, in one of the most problematic 
case of intra-pixel variations (jitter ‘sinl’). As ex¬ 
pected, binning properties depend on the ampli¬ 
tude of white noise relative to systematic noise; 
therefore, for cases with lower SNR it may superfi¬ 
cially appear that systematics are better removed 
in the detrending process. Fig. [18] reports the 
transit parameter retrieved from detrended light- 
curves; in representative cases, we calculated the 
error bars. 

Fig. HH shows the amplitude of discrepancies 
between the detrended light-curves and the theo¬ 
retical model as a function of astrophysical Pois¬ 
son noise, for pixel-ICA and PCD method with 
different polynomial orders: 

1. For the inter-pixel effects the efficiency of the 
two methods is limited by the astrophysical 
Poisson noise level, except when it is smaller 
than the instrument HFPN, in which cases 
pixel-ICA slightly outperforms PCD. 

2. For the intra-pixel effects, residual system¬ 
atics are clearly present in high SNR cases, 
while for lower SNR the limit is the Poisso- 
nian threshold. Second order PCD method 
is far from optimal, while higher order poly¬ 
nomials in some cases can do better than 
pixel-ICA, but they do not always improve 
the results, as already mentioned in Sec. 14.31 
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Fig. 16.- Top panels: (blue) raw light-curves simulated with ctpsf =1 p.u., inter-pixel quantum efficiency 
variations over 9x9 array of pixels, jitter ‘sinl’, and poissonian noise with (left) SNR = 447, and (right) 
SNR = 224; (red) correspondent detrended light-curves with pixel-ICA. Middle panels: the same for raw 
light-curves simulated with (Jpsf =0.2 p.u., intra-pixel quantum efficiency variations over 5x5 array of 
pixels, and jitter ‘sinl’. Bottom panels: the same for light-curves simulated with ctpsf =0.2 p.u., intra-pixel 
quantum efficiency variations over 5x5 array of pixels, and jitter ‘cosl’. All the light-curves are binned over 
10 points, as in previous figures. 
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Fig. 17. - Left panel: Root mean square of residuals for binning over 1 to 10 points, scaled to their non- 
binned values, obtained for simulations with apsF =0.2 p.u., intra-pixel effects over a 5x5 array, jitter ‘sinl’, 
and poissonian noise with SNR = 447. The dashed black line indicates the expected trend for white residuals, 
blue dots are for normalized raw light-curves, and red V are for pixel-ICA detrendend light-curves. Right 
panel: the same, but with SNR = 224. 
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Fig. 18.- Top panel: best estimates of the planet-to-star radii ratio, p = r p /R Sl for some detrended light- 
curves with Poissonian noise at different levels. Error bars are reported for representative cases. Middle 
panel: the same for the orbital semimajor axis in units of the stellar radius, ao = a/R s . Bottom panel: the 
same for the orbital inclination, i. 
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Fig. 19.- Top panel: Root mean square of discrepancies between detrended light-curves with (red circles) 
pixel-ICA, and (green) PCD method, and the theoretical model, obtained for simulations with crpsF =1 
p.u., inter-pixel effects over a 9x9 array, jitter ‘sinl’, and poissonian noise at different levels. Black “+’ 
along dashed line indicate the root mean square of the poissonian signal. Middle panel: The same for 
simulations with crpsF =0.2 p.u., intra-pixel effects over a 5x5 array, jitter ‘sinl’. Green “+’ refer to light- 
curves detrendend with a 4 th order polynomial, and green rhomboids with a 6 th order one. Bottom panel: 
The same for simulations with crpsF =0.2 p.u., intra-pixel effects over a 5x5 array, jitter ‘cosl’. 
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4.5. Inter-pixel effects without noise 

In Sec. 14.11 we state that inter-pixel effects are 
well detrended with pixel-ICA method, based on 
the binning properties of residuals (and consis¬ 
tent results). Given that binning properties can 
only prove that systematics are negligible com¬ 
pared to the actual white noise level, we performed 
a last test for a simulation with inter-pixel ef¬ 
fects ( (Jpsf =1, 9x9 array, jitter sinl) and a re¬ 
duced white noise level, by a factor of 10. Note 
that it is an extreme low value of 0.5 photon 
counts/pixel/data poinld. Fig. l20l shows the raw 
and detrended light-curves for this simulation, and 
the binning properties of their residuals. Time 
structure is very high for the raw light-curve, but 
it is again well detrended by pixel-ICA. We also 
checked that all the retrieved parameters are con¬ 
sistent with the original values within 1 a. 


2 We do not set this noise exactly to 0, because pixel time 
series with many zeroes would negatively affect the ICA 
detrending. 
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Fig. 20. Left panels: (blue) raw light-curve simulated with crpsF =1 p.u., inter-pixel sensitivity variations 
over 9x9 array of pixels, jitter ‘sinl’, and white noise at 0.5 photon counts/pixel/frame; (red) correspondent 
detrended light-curve with pixel-ICA. Right panel: Root mean square of residuals for binning over 1 to 10 
points, scaled to their non-binned values. 
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4.6. Testing a variant of pixel-ICA algo¬ 
rithm 


In sections 14.1114.21 and 14.31 data are detrended 
before fitting the astrophysical model. In the pre¬ 
vious literature, there are two schools of thought: 
(1) to use the out-of-transi t data only to calibrat e 
the systemati cs mode l (e.g. iBeaulieu et al.l ( 201 ill : 
iBallard et al.1 ( 201dh ') . (2) fitting for systemat- 
ics and astroph y sical m odels simultaneou sly 


e. 


Knutson et al. 


Deming et al 



IStevenson et al.l (I2012D 
A major advantage of 
the second approach is that correlations between 
transit and decorrelating parameters can be in¬ 
vestigated. However, there are two arguments in 
favour of the first approach: 


1. the out-of-transit can be used to calibrate 
the instrumental modulations, minimizing 
the risk of confounding them with the as¬ 
trophysical signal. The risk of confusion is 
higher if the timescale of a systematic signal 
is similar to the transit duration; 

2. possible bumps due to occulted star spots 
during the transit would be reduced or can¬ 
celed by an erronoeous correction of the in¬ 
strument systematics. 


The second argument does not apply to our simu¬ 
lations, given that the only astrophysical signal is 
the planetary transit. In general, the best option 
is to check whether the two approaches lead to the 
same results, and if not, to investigate the origin 
of the discrepancy. 

A variant of the original pixel-ICA algorithm 
consists in fitting the independent components 
plus a transit m odel to the raw ligh t-curves (see 
Appendix C.5 in lMorello et alJ ( 2015 1). We tested 


that, in our simulations, this alternative algorithm 
essentially leads to the same results obtained with 
the original algorithm. A few exceptions occur for 
some of the configurations obtained with ‘cosl’ jit¬ 
ter (unfortunate timescale and phasing), for which 
results obtained with the two approaches are not 
identical, but consistent within ler. In all tested 
cases, the (partial) parameter error bars obtained 
by the MCMC chains are very similar for the two 
algorithms, which indicates the absence of strong 
correlations between decorrelating and transit pa¬ 
rameters. However, we suggest that the final er¬ 
ror bars should be computed according to Eq. [Tj 
to account for intrinsic errors on the components 
extracted. 


4.7. Comparison between 
PLD algorithms 

Deming et al 


pixel-ICA and 


Recently, 


( 20141 ) proposed a 
pixel-level detrending method. The underlying 
model to adapt to the data is the product of an 
astrophysical factor, a generalized function of the 
pixel intensities, and a temporal term: 


S f = 5.(1- DE(t)) x F{Pl Pi ..., P 4 ) x Q(t) (6) 


where 5* is the (constant) stellar luminosity, D is 
the transit/eclipse depth, E(t) is the eclipse shape 
normalized to unit amplitude, P\ is the number 
of photon counts read in pixel i at time t. First 
order expansions of the the three factors in Eq. [6] 
lead to the following formula: 


AS 4 = c t SP! + DE(t) + ft + gt 2 + h (7) 

i= 1 


where A indicates the total fluctuations from all 
sources, c,; are partial derivatives from the Tay¬ 
lor expansion, (ft + gt 2 ) represent the temporal 
variatio ns (it is possib l e to a dopt other functional 
forms). iDeming et al.l ( 2014 ) recommend the fol¬ 
lowing replacement: 


5Pl = P* = (8) 

Fitting Eq. [7]to our simulated light-curves, the al¬ 
gorithm did not converge. A possible explanation 
is that the cross-terms between pixel fluctuations 
and transit are not negligible. Therefore, we read- 
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justed Eq. [7] to include those terms: 

*S t =(£c i Pt + hj(l-DE(t)) (9) 

The explicit temporal dependence is ignored, given 
that it is not included in our simulations. 

We fitted for three transit parameters (p, ao 
and i), and Cj and h coefficients simultaneously. 
Results obtained with this method are practically 
indistinguishable from the ones obtained with the 
latest variant of pixel-ICA. 


4. In most cases, pixel-ICA outperforms PCD 
method, especially if the instrument PSF is 
narrow, or it is not entirely within the pho¬ 
tometric aperture. 

5. Intra-pixel effects are only detectable if the 
PSF is relatively small. 

6. Intra-pixel effects cannot be totally de¬ 
trended by any of the three methods, but 
pixel-ICA, in most cases, outperforms PCD 
method, which is largely case-dependent. 
Also, pixel-ICA method provides consistent 
results within the error bars. 


5. Conclusions 


We have tested the pixel-ICA algorithm, i .e. a 


non-p a ramet ric method proposed bv lMorello et al 


(2014, 12015 ) to detrend Spitzer/IRAC primary 
transit observations, on simulated datasets. Sys- 
tematics similar to the ones present in Spitzer/IRAC 
datasets are obtained by combining instrumental 
jitter with inter- or intra-pixel sensitivity varia¬ 
tions. A variety of jitter time series is used to test 
the pixel-ICA method with: 


1. periodic signals with different frequencies, 
phasing, and shape; 

2. non-stationary signals with varying ampli¬ 
tudes or frequencies; 

3. sudden change. 


The detrending performances of pixel-ICA method 
have been compared with division by a polyno¬ 
mial function of the centroi d, in this paper PCD 
method, and PLD method ( Deming et al • Mi). 
Here we summarize the main results found: 


1. Pixel-ICA algorithm can detrend non-stationary 
signals and sudden changes, as well as pe¬ 
riodic signals with different frequencies and 
phasing, relative to the transit. 

2. Inter-pixel effects are well-detrended with 
pixel-ICA method. 

3. Even if the instrument PSF is not entirely 
within the array of pixels, pixel-ICA leads 
to results which are consistent at ~ler with 
the input parameters. 


7. It is possible to fit the astrophysical signal 
after detrending or together with the other 
components. The only differences are regis¬ 
tered if at least one of the non-transit com¬ 
ponents has a similar shape at the time of 
transit, in which case the first approach is 
preferable, but the two results were consis¬ 
tent within lcr. 

8. The PLD method, updated to include cross¬ 
term between pixel fluctuations and the as¬ 
trophysical signals, lead to very similar re¬ 
sults than pixel-ICA, particularly if the as¬ 
trophysical signal is fitted together with the 
other components. 


In conclusion, we have found in a variety of sim¬ 
ulated cases that pixel-ICA performs as well or 
better than other methods used in the litera¬ 
ture, in particular polynomial centroid correc¬ 


ti ons and pixel-level decorrelation (iDemin g et al 


20141) . The main advantage of pixel-ICA over 


other methods relies on its purely statistical foun¬ 
dation without the need of imposing prior knowl¬ 
edge on the instrument systematics, therefore 
avoiding a potential source of error. The results 
of this paper, together with previ ous analyses of 
real Spitzer/IRAC datasets ( Morello et al. 2014 . 
2015h . suggest that photometric precision and re¬ 
peatability at the level of one part in 10 4 can be 
achieved with current infrared space instruments. 


The author would like to thank Prof. G. Tinetti 
and Dr. I. P. Waldmann for useful comments. G. 
Morello is funded by UCL Perren/Impact scholar¬ 
ship (CJ4M/CJ0T). This work was partially sup¬ 
ported by ERC project numbers 617119 (Exo- 
Lights). 
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A. Independent vs Principal Component Analysis 

Principal Component Analysis (PCA) is a statistical technique to separate various components in a set 
of observations/recordings. The main difference between PCA and ICA is that the former just assumes the 
components to be uncorrelated, while the criterion of separation adopted in ICA is the maximum statistical 
independence between the components, which is a more constraining property. 

Two scalar random variables x and y are uncorrelated if their covariance is zero, or equivalently: 

E{xy} = E{x}E{y} (Al) 


where !£{•} denotes the expected value. 

Two variables are statistically independent if their joint probability density factorizes into the product of 
their marginal densities: 

Px, v {x,y) =Px{x)p v (y) (A2) 

Equivalently, they must satisfy the property: 

E{f{x)g{y)} = E{f{x)}E{g(y )} (A3) 

being f{x) and g(x) any integrable functions of x and y, respectively. 

Therefore, independent variables are also uncorrelated, but not viceversa. Only if the variables have gaussian 
distributions, uncorrelatedness implies independence. 

In practice: 


• PCA uses up to second-order statistics. It can only deals with mixtures of gaussian signals, where 
uncorrelatedness equals independence. 

• ICA uses third and fourth order statistics, and it can separate mixtures of non-gaussian signals. 


Many ICA algorithms ran PCA as a useful preprocessing step 


(jHvvarinen et al. 2001). 
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B. Additional tables 
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Table 7: Retrieved transit parameters for simulations with apsF =1, 9x9 array, inter-pixel effects (see Sec. 
ED - In representative cases, we report the partial error bars obtained by the residuals, the final error bars, 
and the worst case error bars (see Sec. [2j . 


Jitter 

Parameters 

Best values 

l-cr errors 

l-cr errors 

l-cr errors 




(residual scatter only) 

(ICA) 

(ICA worst case) 


V 

0.15505 

incur 5 

1.4xl0 -4 

1.4xl0 -4 

sinl 

do 

8.99 

0.02 

0.03 

0.04 


i 

85.78 

0.03 

0.04 

0.04 


V 

0.15501 

l.OxKP 4 

1.5xl0 -4 

1.5xl0 -4 

cosl 

O-a 

9.05 

0.03 

0.04 

0.04 


i 

85.85 

0.03 

0.04 

0.05 


V 

0.15507 




sin2 

Oq 

8.98 





i 

85.78 





V 

0.15510 




cos2 

do 

8.94 





i 

85.74 





V 

0.15505 




sin3 

do 

8.99 





i 

85.79 





V 

0.15503 




cos3 

do 

9.02 





i 

85.82 





V 

0.15506 




sawl 

do 

8.98 





i 

85.78 





V 

0.15506 




sawlvl 

do 

8.99 





i 

85.79 





V 

0.15511 




sawlv2 

do 

9.02 





i 

85.81 





V 

0.15508 

l.OxKT 4 

1.5xl0 -4 

1.5 xlO -4 

sawlv3 

do 

8.99 

0.03 

0.04 

0.04 


i 

85.79 

0.03 

0.05 

0.05 


V 

0.15509 




sawlvfl 

do 

8.98 





i 

85.77 





V 

0.15506 

l.OxKP 4 

1.5xl0 -4 

i.6xur 4 

sawlvf2 

do 

8.99 

0.03 

0.04 

0.05 


i 

85.79 

0.03 

0.05 

0.05 


V 

0.15503 

l.OxlO -4 

1.4xl0 -4 

1.5xl0 -4 

jump04c 

do 

9.00 

0.03 

0.04 

0.04 


i 

85.80 

0.03 

0.04 

0.05 
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Table 8: Retrieved transit parameters for simulations with <jpsf =1, 5x5 array, inter-pixel effects (see Sec. 
ED - In representative cases, we report the partial error bars obtained by the residuals, the final error bars, 
and the worst case error bars (see Sec. [2j . 


Jitter 

Parameters 

Best values 

l-o errors 

l-o errors 

1-cr errors 




(residual scatter only) 

(ICA) 

(ICA worst case) 


V 

0.15524 

l.sxicr 4 

2.2x10-“ 

2.3x10-“ 

sinl 

do 

8.97 

0.04 

0.06 

0.07 


i 

85.76 

0.05 

0.07 

0.07 


V 

0.15505 

1.6xl0~ 4 

2.2x10-“ 

2.3x10-“ 

cosl 

O-a 

9.18 

0.05 

0.07 

0.07 


i 

85.99 

0.05 

0.07 

0.07 


V 

0.15510 




sin2 

Oq 

9.08 





i 

85.88 





V 

0.15536 




cos2 

do 

8.89 





i 

85.66 





V 

0.15528 




sin3 

do 

8.99 





i 

85.79 





V 

0.15516 




cos3 

do 

9.07 





i 

85.86 





V 

0.15519 




sawl 

do 

8.97 





i 

85.77 





V 

0.15517 




sawlvl 

do 

9.02 





i 

85.81 





V 

0.15494 




sawlv2 

do 

8.92 





i 

85.74 





V 

0.15527 

l.sxiir 4 

2.4x10-“ 

2.6x10-“ 

sawlv3 

do 

8.98 

0.04 

0.07 

0.07 


i 

85.77 

0.05 

0.07 

0.08 


V 

0.15524 




sawlvfl 

do 

9.00 





i 

85.78 





V 

0.15531 

i.6xur 4 

2.2x10““ 

2.4x10-“ 

sawlvf2 

do 

8.99 

0.04 

0.06 

0.07 


i 

85.79 

0.05 

0.07 

0.07 


V 

0.15508 

1.3xl0 -4 

2.6x10-“ 

2.8x10-“ 

jump04c 

do 

9.04 

0.04 

0.07 

0.08 


i 

85.83 

0.04 

0.08 

0.09 
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Table 9: Retrieved transit parameters for simulations with upsF =0.2, 5x5 array, inter-pixel effects (see Sec. 
14.21) . In representative cases, we report the partial error bars obtained by the residuals, the final error bars, 
and the worst case error bars (see Sec. [2j . 


Jitter 

Parameters 

Best values 

1-a errors 

l-<7 errors 

1-cr errors 




(residual scatter only) 

(ICA) 

(ICA worst case) 


V 

0.15507 

3x10-* 

6x10-° 

7xl0~“ 

sinl 

clq 

9.010 

0.010 

0.018 

0.020 


i 

85.808 

0.010 

0.020 

0.022 


V 

0.15505 

4xl0 _!> 

6x10-“ 

7xl0~“ 

cosl 

clq 

9.018 

0.010 

0.016 

0.019 


i 

85.815 

0.011 

0.017 

0.020 


V 

0.15508 




sin2 

CLQ 

9.009 





i 

85.806 





V 

0.15505 




cos2 

do 

9.004 





i 

85.801 





V 

0.15510 




sin3 

CLQ 

9.001 





i 

85.798 





V 

0.15507 




cos3 

do 

9.015 





i 

85.812 





V 

0.15506 




sawl 

CLQ 

8.993 





i 

85.791 





V 

0.15506 




sawlvl 

do 

9.000 





i 

85.798 





V 

0.15506 




sawlv2 

do 

8.997 





i 

85.795 





V 

0.15506 

3XKP 5 

5xl0-“ 

6xl0~“ 

sawlv3 

do 

8.998 

0.009 

0.015 

0.018 


i 

85.796 

0.010 

0.016 

0.019 


V 

0.15508 




sawlvfl 

do 

8.999 





i 

85.797 





V 

0.15506 

3X10- 5 

5x10““ 

6x10““ 

sawlvf2 

do 

9.001 

0.009 

0.015 

0.018 


i 

85.798 

0.010 

0.016 

0.019 


V 

0.15508 

— 

X 

o 

6x10-“ 

6xl0~“ 

jump04c 

do 

9.00 

0.011 

0.017 

0.019 


i 

85.80 

0.011 

0.018 

0.020 
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Table 10: 
Sec. 14.31) . 
bars, and 


Jitter 

Parameters 

Best values 

l-cr errors 

l-<7 errors 

l-cr errors 




(residual scatter only) 

(ICA) 

(ICA worst case) 


V 

0.1551 

3xl0 -4 

4xl0 -4 

5xl0“ 4 

sinl 

do 

9.05 

0.07 

0.10 

0.13 


i 

85.84 

0.08 

0.11 

0.14 


V 

0.1550 

2xl0 -4 

3xiir 4 

3xl0 -4 

cosl 

O-a 

9.12 

0.07 

0.08 

0.09 


i 

85.92 

0.08 

0.08 

0.09 


V 

0.1550 




sin2 

Oq 

9.13 





i 

85.93 





V 

0.1551 




cos2 

do 

9.10 





i 

85.89 





V 

0.1551 




sin3 

do 

9.01 





i 

85.81 





V 

0.1551 




cos3 

do 

9.02 





i 

85.82 





V 

0.1548 




sawl 

do 

9.10 





i 

85.93 





V 

0.1555 




sawlvl 

do 

8.92 





i 

85.71 





V 

0.1546 




sawlv2 

do 

9.27 





i 

86.11 





V 

0.1555 

4xl0 -4 

4xlCr 4 

5xl0“ 4 

sawlv3 

do 

8.95 

0.10 

0.11 

0.13 


i 

85.77 

0.11 

0.12 

0.15 


V 

0.1555 




sawlvfl 

do 

8.78 





i 

85.54 





V 

0.1550 

3xl0 -4 

3xl0 -4 

3xl0 -4 

sawlvf2 

do 

8.99 

0.09 

0.10 

0.10 


i 

85.79 

0.10 

0.11 

0.11 


V 

0.15508 

— 

X 

o 

3xl0 -4 

4xl0 -4 

jump04c 

do 

9.001 

0.011 

0.10 

0.12 


i 

85.796 

0.012 

0.10 

0.13 


Retrieved transit parameters for simulations with <jpsf =0.2, 5x5 array, intra-pixel effects (see 
In representative cases, we report the partial error bars obtained by the residuals, the final error 
the worst case error bars (see Sec. [2]). 
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C. List of acronyms 

HFPN = High-Frequency Pixel Noise 
ICA = Independent Component Analysis 
ISR = Interference-to-Signal Ratio 
MCMC = Markov Chain MonteCarlo 
PCA = Principal Component Analysis 
PCD = Polynomial Centroid Division 
PLD = Pixel-Level Decorrelation 
PSF = Point Spread Function 
RMS = Root Mean Square 
SNR = Signal-to-Noise Ratio 
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